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ABSTRACT 

Properties, structure, and thermal evolution of neutron stars are determined 
by the equation of state of stellar matter. Recent data on isospin-diffusion and 
isoscaling in heavy-ion collisions at intermediate energies as well as the size of 
neutron skin in ^osp^ have constrained considerably the density dependence of 
the nuclear symmetry energy and, in turn, the equation of state of neutron-rich 
nucleonic matter. These constraints could provide useful information about the 
global properties of rapidly rotating neutron stars. Models of rapidly rotating 
neutron stars are constructed applying several nucleonic equations of state. Par- 
ticular emphasis is placed on configurations rotating rigidly at 716 and 1122Hz. 
The range of allowed hydrostatic equilibrium solutions is determined and tested 
for stability. The effect of rotation on the internal composition and thermal prop- 
erties of neutron stars is also examined. At a given rotational frequency, each 
equation of state yields a range of possible neutron stars configurations restricted 
by the Keplerian (mass-shedding) limit, corresponding to the maximal circum- 
ferential radius, and the limit due to the onset of instabilities with respect to 
axial-symmetric perturbations, corresponding to the minimal equatorial radius 
of a stable neutron star model. We show that the mass of a neutron star rotating 
uniformly at 1122Hz is between 1.7 and 2.1Mq. Central stellar density and pro- 
ton fraction decrease with increasing rotational frequency with respect to static 
models, and depending on the exact stellar mass and angular velocity, can drop 
below the Direct Urea threshold thus closing the fast cooling channel. 



Subject headings: dense matter — equation of state — stars: neutron — stars: 
rotation 
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Introduction 



Neutron stars are one of the most exotic objects in the universe. Matter in their cores 
is compressed to huge densities ranging fro m the density of no rmal nuclear matter, po 
0.16/m~^, to an order of magnitude higher (iGlendenningj |2000| ) . The number of baryons 
forming a neutron star is in the order of A ^ 10^^. Understanding properties of matter under 
such extreme conditions of density (and pressure) is still far from complete and represents 
one of the most important but also challenging problems in modern physics. 

S trictly speak ing, neutron stars are associated with two classes of astrophysical ob- 
jects ( IWeberlll999l ). Pulsars belong to the first class and are generally accepted to be rotating 
neutron stars. X-ray sources belong to the second class and some of them are neutron stars 
on close binary orbit with an ordinary star (e.g.. Her X-1 and Vela X-1) (IWeberl 119991). Be- 



cause of their strong gravitational binding neutron stars can rotate very fast (IBeiger et al. 
2007 ). The first millisecond pulsar PSR1937-I-214, spinning at z/ = QAlHz (IBacker et al. 



19821 ). was discovered in 1982, and during the next decade or so almost every year a new 



one was reported. In the recent years the situation changed considerably with the discover y 
of an anomalously large population of millisecond pulsars in globular clusters (IWeberlll999l ). 
where the density of stars is roughly 1000 times that in the field of the galaxy and which are 
therefore very favorable sites for formation of rapidly rotating neutron stars which have been 
spun up by the means of mass accretion from a binary companion. Presently the number of 
the observed pulsars is close to 2000, and the detection rate is rather high. 



In 2006 lHessels et al.l (120061 ) reported the discovery of a very rapid pulsar J1748-2446ad, 
rotating a.t u = 71QHz and thus breaking the previous record (of QAlHz). However, even 
this high rot ational frequency i s too low to affect the structure of neutron stars with masses 
above IMq (iBejger et al.ll2007l ). Such pulsars belong to the slow-rotation regime since their 
frequencies are considerably lower than the Kepler (mass-shedding) frequency u^. (The mass- 
shedding, or Kepler, frequency is the highest possible frequency for a star before it starts to 
shed mass at the equator.) Neutron stars with masses above IMp) enter the rap id-rotation 
regime if their rotational freq uencies are higher than lOOOHz (jBejger et al.l 120071 ). A recent 
report by lKaaret et al.l (120071 ) suggests that the X-ray transient XTE J1739-285 contains the 
most rapid pulsar ever detected rotating at z/ = 1122Hz. T his discovery has re awaken the 
interest in building models of rapidly rotating neutron stars (jBejger et al.l 120071 ). 



Since neutron stars are objects of extremely condensed matter, the geometry of space- 
time is considerably altered from that of a flat space. Therefore, the construction of 
realistic models of neutron stars has to be done in the framework of General Relativ- 



ity (IWeberl 1 19991 : IGlendenningj |2000| ). Detailed knowledge of the equation of state (EOS) 
of stellar matter over a very wide range of densities is required for solving the neutron- 
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star structure equations. At present time the behavior of matter under extreme densi- 
ties such as those found in the interiors of neutron stars is still highly uncertain and re- 
lies upon, often, rather controversial theoretical predictions. Fortunately, heavy-ion reac- 
tions provide a unique means to con strain the EOS of dense nuclear matter in terrestrial 
laboratories (jPanielewicz et al.l |2002| ). In particular, one of the main sources of uncer- 
tainty in__tlie_jtelki;_jnB^ en- 



ergy (ILattimer fc Prakashll2004l : ISteiner et al.ll2005l : iKrastev fc Sammarrucall2006l ). esym{p) 



which is the difference between the energies of pure neutron and symmetric nuclear matter. 
Due to its importance for the neutron star structure determining the density dependence 
of the nuclear symmetry energy has been a priority goal for the intermediate energy heavy- 
ion community. Although extracting the symmetry energy is not an easy task due to the 
complicated role of the isospin degree of freedom in reaction 
prob es of the sy mmetry energy have been suggested (Li et al. 19971: 



19981 ) (see also (ILi fc Udo Schroederl 1200 ll : iDanielewicz et al. 



2002 



Li 


2000|. 


2002 


■ 


Li et al. 


Bar an et al. 


20051) for 



reviews). Some significant progress have been made recently in determining the density 
depen dence of Cgym a t subsaturat i on densities using: (1) isospin diffusion (ITsang et al. 



20041 ) and isoscaling (ITsang et al. 



react i ons at intermedia te energies (IShi fc Danielewicz 



2001 



Shettv. Yennello fc Souliotis 



2003 



20071) in heavy-ion 



Chen et al. 



20051: 



Li &: Chen! 120051) . and ( 2) sizes of neutron skins in heavy nuc lei 



200,4 



Steiner 



Steiner et al. 



& Li 



20051 : 



Horowitz &: Piekarewical200ll . 120021 : iTodd-Rutel fc Piekarewiczll2005l ). At supranormal den- 
sities, a number of potential probes of the symnietry energy have been proposed although 
there is not much data available yet (jChen et al.l 120071 ). 
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2003) for a review). In this paper we com 





bine recently obtained data on isospin diffusion, information from flow observables, studies 
of neutron skin data of ^ospft and other information to constrain global properties of rotating 
neutron stars. Applying several nucleonic equations of state and the i?A^S0 code developed 



-'^ Thanks to Nikolaos Stergioulas the RNS code is available as a public domain program at 
http:/ /www-gravity, phys.uwm.edu/rns/i 
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and made available to the public by Nikolaos Stergioulas (IStergioulas fc FriedmanI Il995l ) , 
we construct one-parameter 2-D stationary configurations of rapidly rotating neutron stars. 
The computation solves the hydrostatic and Einstein field equations for mass distributions 
rotating rigidly under the assumptions of stationary and axial symmetry about the rotational 
axis, and reflection symmetry about the equatorial plane. This work is organized in the fol- 
lowing way: after the introduction in Section 1, we review briefly the structure equations 
of both static and rotating neutron stars on Section 2. Our numerical results are presented 
and discussed in Section 3. We conclude with a short summary in Section 4. 



2. Equations of neutron star structure 



In what follows we present a brief review of structure equations of both static and 
rotating neutron stars. As already mentioned in the introduction neutron stars are objects of 
extremely compressed matter and therefore proper understanding of their properties requires 
application of both General Relativity and the theories of dense matter, which constitute 
nuclear and particle physics problem. In this respect neutron stars provide a direct link 
between two of t he froiitiers of modern physics - General Relativity and strong interactions 
in dense matter (jWeberlll999l ). The connection between both branches of physics is provided 
by Einstein's field equations 



Qt^y = ^t^y _ = 87rT^'^(e, P(e)), 



1^ 



(/i, 1/ = 0, 1, 2, 3) which couple the Einstein curvature tensor, G^^ , to the energy-momentum 
tensor, 

T^^ = (e + P)u^'u'' + P^'^^ (2) 

of stellar matter. In the above equations P and e denote pressure and mass energy density, 
while R^^ , g'^", and R denote the Ricci tensor, the metric tensor, and the Ricci scalar 



curvature respectively (see e.g.. iGlendenningl 120001 ). In Eq. ([2]) is the unit time-like four- 
velocity satisfying u'^u^ = —1. The tensor T^'^ contains the EOS of stellar matter in the 
form P(e). In general, Einstein's field equations and those of the nuclear many-body problem 
were to be solved simultaneously since the baryons and quarks follow the geodesies of the 
curved space-time whose geometry, determi ned by the Ei nstein's field equations, is coupled 
to the total mass energy density of matter (jWeberlll999l ). In the case of neutron stars, as 
for all astrophysical situations for which the long-range gravitational force can be separated 
from the short-range strong force, the deviation from flat space-time at the length-scale of 
the strong interactions (~ is practically zero up to the highest densities achieved in 

the neutron star interiors. (This is not to be confused with the global length-scale of neutron 
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stars (~ 10km) for which M/R ~ 0.3 depending on the star's mass (in units c = G = 1 
so that Mq 1.475A;m).) In other words, gravity curves space-time only on a macroscopic 
scale but to a very good approximation leaves it flat on a microscopic scale. To achieve 
an appreciable curvature on a microscopic level at which the strong interactions do minate 



the particle dynarn ics mass densities greater than ~ 10 g cm would be necessary (1 Weber 



19991 : lThornelll966l ) . Under this circumstances the problem of constructing models of neutron 
stars separates into two distinct tasks. First, the short-range effects of the nuclear forces 
are described by the principles of many-body nuclear physics in a local inertial frame (co- 
moving proper reference frame) in which space-time is fiat. Second, the coupling between the 
long-range gravitational force and matter is accounted for by solving the general relativistic 
equations for the gravitational field described by the curvature of space-time, leading to the 
global structure of stellar configurations. 



2.1. Static stars 

In the case of spherically symmetric static (non-rotating) stars the metric has the famous 
Schwarzschild form: 



ds^ = -e^^i-Ut' + e^^^'-^dr^ + r\d9' + sin^ 9dc^') 
c = G = 1) where the metric functions (f){r) and A(r) are given by: 

e2AW = (l-7(r))-\ 

1 - 7(r)) r > Rstar, 

2m(r) 



g2<^(r) ^ g-2A(r) 



with 



ii r < R 



star 



7(r) 



(3) 

(4) 
(5) 

(6) 



^ if r > i? 



■star 



For a static star Einstein's field equat ions (Eg. ([J)) reduce then to the farniliar Tolman- 
Oppenheimer-Volkoff equation (TOV) (jTolmarull939l : lOppenheimer fc Volkofllll939l ): 



dP{r) e(r)m(r) 
dr 



1 



P(r) 



e r 



1 H 9-^ 



1 - 



2m(r) 



-1 



(7) 



where the gravitational mass within a sphere of radius r is determined by 

dm{r) 



dr 



47re(r)r dr 
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The metric function 0(r) is determined through the following differential equation: 

d(j){r) _ 1 dP{r) 

dr e(r) + P(r) dr ' 

with the boundary condition at r = i? 



(9) 



0(r = /2) = ^ln(l-^) (10) 



To proceed to the solution of these equations, it is necessary to provide the EOS of 
stellar matter in the form P{e). Starting from some central energy density ec = e(0) at the 
center of the star (r = 0), and with the initial condition m(0) = 0, the above equations 
can be integrated outward until the pressure vanishes, signifying that the stellar edge is 
reached. Some care should be taken at r = since, as seen above, the TOV equation is 
singular there. The point r = R where the pressure vanishes defines the radius of the star 
and M = m{R) = An e{r')r''^dr' its gravitational mass. 

For a given EOS, there is a unique relationship between the stellar mass and the central 
density ec- Thus, for a particular EOS, there is a unique sequence of stars parameterized by 
the central density (or equivalently the central pressure -P(O)). 



2.2. Rotating stars 



Equations of stellar structure of (rapidly) rota ting neutron stars are considerably more 
complex than those of spherically symmetric stars (IWeberlll999l ). These complications arise 
due to the rotational deformations in rotating stars (i.e., flattening at the poles and bulging 
at the equator), which lead to a dependence of the star's metric on the polar coordinate 6. 
In addition, rotation stabilizes the star against gravitational collapse and therefore rotating 
neuron stars are more massive than static ones. A larger mass, however, causes greater 
curvature of space-time. This renders the metric functions frequency-dependent. Finally, 
the general relativistic effect of dragging the local inertial frames implies the occurrence of 
an additional non-diagonal term, g'^'f', in the metric tensor g'^'^. This term imposes a self- 
consistency condition on the stellar structure equations, since the degree at which the local 
inertial frames are dragged along by the star, i s determined by the initially unknown stellar 
properties like mass and rotational frequency (IWeberlll999l ). 



Here we outline briefly the equations solved by the RNS code. The coordinates of the 
stationary, axial symmetric space-time used to model a (rapidly) rotating neutron star are 



- 7- 



defined through a generahzation of Bardeen's metric (IStergioulaa l2003l ) : 

+ e"'-Pr^sm^e{d(f)-ujdt^), 



where the metric potentials 7, p, a, and the angular velocity of the stellar fluid relative to 
the local inertial frame, u, are functions of the quasi-isotropic radial coordinates, r, and the 
polar angle 6 only. Th e matter inside a rigidly rotating star is approximated as a perfect 
fluid (jStergioulasll2003l ) . whose energy momentum tensor is given by Eq. ([2]). The pro per ve- 
locity of matter, v, relative to the local Zero Angular Momentum Observer (ZAMO) (lOuyed 
2002h is defined as 

(12) 



V 



uo e 



-p(r) 



with 



rsin(^)(r] 

/vP the angular velocity of a fluid element. The four- velocity is given by 



-(7+P)/2 



;i,o,o,fi) 



(13) 



In the above equation the function (7 + p)/2 represents the relativistic generalization of the 
Newtonian gravitational potential, while exp[(7 + p)/2] is a time dilation factor between 
an observer moving with angular velocity uj and one at infinity. Substitution of Eq. f[T^ 
into Einstein's fields equations projected onto the ZAMO reference frame gives three elliptic 
partial differential equations for the metric potentials 7, p, and tu, and two linear ordinary 
differential equations for the metric potential a. Technically, the elliptic differential equations 
for the metric functions are converted into integral equations which are then s olved iteratively 
applying Green's function approach (IKomatsu et al.lll989l : IStergioulasll2003l ). 



From the relativistic equations of motion, the equations of hydr ostatic equilibrium for 
a barotropic fluid may be obtained as (IStergioulad l2003l : IOuyedll2002l ): 



h{P) -hp = -[ujp + pp--/ - p- ln(l - v^) + A\uj - n,) 



(14) 



with h{P) the specific enthalpy, Pp the re-scaled pressure, hp the specific enthalpy at the 
pole, 7p a nd Pp the values o f the metric potentials at the pole, flc = ''"e^, and A a rotational 
constant (IStergioulasI |2003| ) . The subscripts p, e, and c label the corresponding quantities 
at the pole, equator and center respectively. The RNS code solves iteratively the integral 
equations for p, 7 and uj, and the ordinary differential equation for the metric function a 
coupled with Eq. (IMI) and the equations for hydrostatic equilibrium at the stellar center and 
equator (given h{Pc) and h{Pe) = 0) to obtain p, 7, a, u, the equatorial coordinate radius 
Te, angular velocity fl, energy density e, and pressure P throughout the star. 



- 8 - 




Fig. 1. — (Color online) Equation of state. The upper frame shows the mass energy density 
as a function of baryon density (in units of po) and the lower frame shows the total pressure 
(including the lepton contributions) versus baryon density. (The "dip" exhibited by the 
density curve of the APR EOS is due to a phase tra nsition from low density phase (LDP) 
to high density phase (HDP). See f Akmal et al.lllQQsl ) for details.) 



3. Results and discussion 



We compute properties of (rapidly) rotating neutron stars employing several nucleonic 
EOSs and the RNS code. The EOSs applied here are shown i n Fig. 1. We p ay partic- 
ular attention to the EOS computed with the MDI interaction (IDas et al.l l2003l ) since its 
symmetry energy is constrained in the sub-saturation density region by the available nu- 
clear laboratory data. The EOS of symmetric nuclear matter for the MDI interaction is 
constrained up to five times the normal nuclear matter density by the available data on 
collective flow in relativistic heavy-ion collisions. The parameter x is introduced in the 
single-particle potential of the MDI EOS to reflect the largely uncertain density dependence 
of the nuclear symmetry e nergy e,^ymi p ) as p redic ted by various man y-body approaches. 
Since, as demonstrated by iLi fc ChenI (120051 ) and iLi &: Steinerl (120061 ) . only equations of 
state with x between -1 and have symmetry energies consistent with the isospin diffu- 
sion data and measurements of the skin thickness of ^°^Pfe, we thus consider only these 
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two limiting cases as boundaries of the possible rotating neutron star models. Moreover, 
it is interesting to note that the symmetry energy extracted very recently from the isoscal- 
i ng analyses of heavy-ion reaction s is consistent with the MDI calculation using x = 



(IShetty. Yennello &: Sou 



the neutron-star radius ( 



iotis 



20071). The MDI EOS has been recently applied to constrain 



Li fc Steineiil2006l ) with a suggested range compatible with the best 



estimates from observations. In addition, it has been also use d to constrain a possible time 
variation of the gravitation al constant ( ? (IKrastev fc Lil 120071 ) via the gravitoc hemical heat 



i ng for malism developed by lJofre et al.l (120061 ). In Fig. 1 we also show results by lAkmal et al. 



( 119981 ) with the A18 + 6v + UIX* interaction (APR) and recent Dirac-Brueckner-Ha rtree- 
Fock (DBHF) calculations ( Alonso fc Sammarruca 20031 : Krastev &: Sammarruca 2006) wit h 
Bonn B One-Boson- Exchange (OBE) potential (DBHF-fBonn B) jMachleidt et al.lll987[ ). 
Below the baryon density of approximately 0.07/m~^ the equations of state shown in Fig. 1 
are supplemented by a crust al EOS, which is more suitab l e for the low density regime. 
Namely, w e apply the EOS by [Pethick. Ravenhall fc Loren2 ( 1995 ) for the inner crust and 
the one by lHaensel &: PichonI (Il994j ) for the outer crust. At the hig hest densities we assume 
a con tinuous functional for the EOSs employed in this work. (See (IKrastev fc Sammarruca 



20061 ) for a detailed description of the extrapolation procedure for the DBHF-|-Bonn B EOS.^ 



The saturation properties of the nuclear equations of state used in this paper are summarized 
in Table 1. 



Table 1: Saturation properties of the nuclear EOSs (for symmetric nuclear matter) shown in 
Fig. 1. 



EOS 




EsiMeV) 


K{MeV) 


esym{po){MeV) 


m*{pQ)/m 


MDI(x=0) 


0.160 


-16.08 


211.00 


31.62 


0.67 


MDI(x=-l) 


0.160 


-16.08 


211.00 


31.62 


0.67 


APR 


0.160 


-16.00 


266.00 


32.60 


0.70 


DBHF+Bonn B 


0.185 


-16.14 


259.04 


33.71 


0.65 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities at the nuclear saturation density: saturation (baryon) density; energy-per-particle; com- 
pression modulus; symmetry energy; nucleon effective mass to average nucleon mass ratio (with 
m = 938.926MeF c"^). 
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Fig. 2. — (Color online) Neutron star masses and radii. Neutron star equatorial radii (upper 
panel) and total gravitational mass (lower panel) versus central energy density ec- Both 
static (solid lines) and Keplerian (broken lines) models are shown. 

3.1. Keplerian (and static) sequences 

In what follows we examine the effect of ultra-fast rotation at the Kepler frequency on the 
neutron star gravitational masses and radii. The equilibrium configurations for both static 
and (rapidly) rotating neutron stars are parameterized in terms of the central mass energy 
density, ec = e(0) (or equivalently central pressure, Pc = -P(O)). This functional dependence 
is shown in Fig. 2, where we display the stellar equatorial radius (upper frame) and total 
gravitational mass (lower frame) versus central energy density for the EOSs applied in this 
work. Predictions for both static and maximally rotating models are shown. The main 
feature of the mass-density plot is that there exists a maximum value of the gravitational 



mass of a neutron star that a given EOS can support (see e.g., IWeberl Il999l ). This holds 
for both static and (rapidly) rotating stars. The sequences shown in Fig. 2 terminate at 
the "maximum mass" point. Comparing the results for static and rotating stars, it is seen 
clearly that the rapid rotation increases noticeably the mass that can be supported against 
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Fig. 3. — (Color online) Mass-radius relation. Both static (solid lines) and Keplerian (broken 
lines) sequences are shown. T he 1 -- tx err or bar corresponds to the measurement of the mass 
and radius of EXO 0748-676 JOzelll2006h . 




1 2 3 

e, (x 10'^gcnT^) 



Fig. 4. — (Color online) Neutron-star polar radius versus central energy density. Both static 
(solid lines) and Keplerian (broken lines) sequences are shown. The labeling of the curves is 
the same as in Fig. 2. 



collapse while lowering the central density of the maximum-mass configuration. This is 
what one should expect, since, as already mentioned, rotation stabilizes the star against 
the gravitational pull providing an extra (centrifugal) repulsion. The rotational effect on the 
mass-radius relation is illustrated in Fig. 3 where the gravitational mass is given as a function 
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of the circumferential radius. For rapid rotation at the Kepler frequency, a mass increase 
up to ~ 17% (Table 3) is obtained, depending on the EOS. The equatorial radius increases 
by several kilometers, while the polar radius decreases by several kilometers (see Fig. 4) 
leading to an overall oblate shape of the rotating star. Table 2 summarizes the properties 
(masses, radii and central energy densities) of the maximum-mass nonrotating neutron star 
configurations. Our studies on the effect of rapid rotation on the upper mass limits for the 
four EOSs considered in the present paper are presented in Table 3. In each case the upper 
mass limit is attained for a model at the mass-shedding limit where u = u^, with central 
density ~ 15% below th at of the static model with the largest mass. These find ings are 
consistent with those by iFriedman et al.l (Il984l ) and IStergioulas fc Friedman! (119951 ) . Table 
3 also provides an estimate of the upper limiting rotation rate of a neutron star. In general, 
softer EOSs permit larger rot ational frequencies sin ce the resulting stellar models are more 
centrally condensed (see e.g.. iFriedman et al.lll984j ). In the last column of Table 3 we show 
the Kepler frequencies computed via the empirical relation 

1/2 ^ D X -3/2 

0.72 ( ^ 1 I 1 (15) 



10%- 



MqJ \l0kmj 



propo sed by IFriedman et al.l (1l989l ). The uncertainty of Eq. (fTSll is ~ 10% (see lHaensel &: Zdunik 
(119891 ) for an i mproved version of th e empirical formula). At the time of constructing the 
above relation Friedman et al.| (1989) did not consider a then unknown class of minimum 
period EOSs (IStergioulaslll996l ) which explains why the numbers in the last column of Table 
3 exhibit larger deviation from the exact numerical solutions (in column five). This is partic- 
ularly pronounced for the APR EOS for which the approximated Kepler frequency deviates 
~ 14% from the exact solution. (Note that the APR EOS has the lowest period among the 
EOSs considered here.) 



Fig. 5 displays the neutron star gravitational mass as a function of the Kepler frequency. 



Table 2: Maximum-mass static (nonrotating) models. 

EOS MmaxiMo) R{km) e,( x lO^^^? cm'^) 



MDI(x=0) 


1.91 


9.89 


3.02 


APR 


2.19 


9.98 


2.73 


MDI(x=-l) 


1.97 


10.85 


2.57 


DBHF+Bonn B 


2.24 


10.88 


2.36 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities for the static models with maximum gravitational mass: gravitational mass; radius; 
central mass energy density. 
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500 1000 1500 2000 
Vk (Hz) 

Fig. 5. — (Color online) Mass versus Keplerian (mass-shedding) frequency z/^. 



z/fc- The vertical doted li nes denote the frequencies of the fastest neutro n stars - XTE J1739- 
285 kaaret et all boOTh (red) and J1748-2446ad Jflessels et all bood ) (green). From this 
figure one can conclude that only neutron stars with masses above approximately 1.7 Mq 
can rotate at z/ = 1122Hz (and/or faster). We discuss this further in the next subsection 
(3.2). 



Table 3: Maximum-mass rapidly rotating models at the Kepler frequency z/ = z/fe. 



EOS 






Increase (%) 


ec(x 10^^(7 cm ^) 






MDI(x= 


=0) 


2.25 


15 


2.59 


1742 


1610 


APR 




2.61 


17 


2.53 


1963 


1699 


MDI(x= 


=-1) 


2.30 


14 


2.21 


1512 


1423 


DBHF^ 


-Bonn B 


2.69 


17 


2.06 


1685 


1510 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities for the maximally rotating models with maximum gravitational mass: gravitational 
mass; its percentage increase over the maximum gravitational mass of static models; central mass 
energy density; maximum rotational frequency; Kepler frequency as computed via the empirical 
relation by. 
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3.2. Rotation at 716 and 1122 Hz 



In this subsection w e study neutron stars rotating at 716 (IHessels et al.ll2006l ) and 1122 
Hz f lKaaret et al.l 120071 ) which are the rotational frequencies of the fastest pulsars as of 
today. Stabihty with respect to the mass-shedding from equator imphes that at a given 
gravitational mass the equatorial ra dius Reg should be smaller than -R^°^ corresponding 
to the Keplerian limit (IBejger et al.ll2007l ). The value of -R^"^ results from the condition 
that the frequency of a test particle at circular equatorial orbit of radius K^"'^ just above the 
equ ator of the actual r otating star is equal to the rotational frequency of the star. As reported 
by (iBejger et al.ll2007l ) the relation between M and Reg at the "mass-shedding point" is very 
well approximated by the expression for the orbital frequency for a test particle orbiting at 
r = Reg in the Schwarzschild space-time created by a spherical mass. The formula satisfying 



Schw. 
orb 



represented by the dotted line in Figs. 6 and 7, is given by 



1 

2^ 



GM 



(16) 



where u = 716Hz in Fig. 6 and u = 1122Hz in Fig. 7 respectively. This formula for the 
Sch warzschild metric c oincides with the one obtained in Newtonian gravity for a point mass 
M JSeiger et al.ll2007h . Eq. ([16]) implies 



Rv 



IAMq 




Fig. 6. — (Color online) Gravitational mass versus circumferential radius for neutron stars 
rotating at u = 716Hz. See text for details. 
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Fig. 7. — (Color online) Gravitational mass versus circumferential radius for neutron stars 
rotating at = 1122Hz. 

with X = 20.94 for rotational frequency u = llQHz (Fig. 6) and x = 15.52 for v = 1122Hz 
(Fig. 7). 

In Figs. 6 and 7 we observe that the range of the allowed masses supported by a given 
EOS for rapidly rotating neutron stars becomes narrower than the one of static configura- 
tions. This effect becomes stronger with increasing frequency and depends upon the EOS. 
For instance, for models rotating at 1122Hz (Fig. 7) for the x=-l EOS the allowed mass 
range is ~ O.IMq. Since predictions from the x=0 and x=-l EOSs represents the limits of 
the neutron star models consistent with the nuclear data from terrestrial laboratories, we 
conclude that the mass of the neutron star in XTE J1739-285 is between 1.7 and 2.1M0. 

3.3. Rotation and proton fraction 

Finally, we study the effect of (fast) rotation on the proton fraction in the star core. In 
Fig. 8 we show the central baryon density (upper frame) and central proton fraction (lower 
frame) as a function of the rotational frequency for fixed-mass models. Predictions from 
both a; = and x = —1 EOSs are shown. We observe that central density decreases with 
increasing frequency. This reduction is more pronounced in heavier neutron stars. Most 
importantly, we also observe decrease in the proton fraction in the star's core. We recall 
that large proton fraction (above ~ 0.14 for rape/i-stars) leads to fast cooling of neutron 
stars through direct Urea reactions. Our results demonstrate that depending on the stellar 




- 16 - 



mass and rotational frequency, the central proton fraction could, in principle, drop below the 
threshold for the direct nucleonic Urea channel and thus making the fast cooling in rotating 
neutron stars impossible. The masses of the models shown in Fig. 8 are chosen so that the 
proton fraction in stellar core is just above the direct Urea limit for the static configurations, 
see Fig. 9 upper and lower frames. The stellar sequences in Fig. 8 are terminated at the 
Kepler (mass- shedding) frequency. In both cases the central proton frac tion drops below th e 
direct Urea limit at frequencies lower than that of PSR J1748-244ad (IHessels et al.ll2006l ). 
This implies that the fast cooling can be effectively blocked in millisecond pulsars depending 
on the exact mass and spin rate. It might also explain why heavy neutron stars (could) 
exhibit slow instead of fast cooling. For instance, with the x = EOS (with softer symmetry 
energy) for a neutron star of mass approximately 1.9Mq, the Direct Urea channel closes at 
u 470Hz. On the other hand, with the x=-l EOS (with stiffer symmetry energy) the 
direct Urea channel can close only for low mass neutron stars, in fact only for masses well 
below the canonical mass of 1.4M0. This is due to the much stiffer symmetry energy (see 
Fig. 9 middle frame) because of which the direct Urea threshold (Fig. 9 lower frame) is 
reached at much lower densities and stellar masses (Fig. 9 upper frame). 
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Fig. 8. — (Color online) Density (upper panel) and proton fraction (lower panel) versus 
rotational frequency for fixed neutron star mass. 
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Fig. 9. — (Color online) Neutron star mass (upper panel), nuclear symmetry energy (middle 
panel) and proton fraction (lower panel) versus baryon density. Predictions from the MDI 
1,0) are shown. 



Before closing the discussion in this section a few comments are in order. In the present 
study we do not consider "exotic" states of matter in neutron stars. On the other hand, 
due to the rapid rise of the baryon chemical potentials several other species of particles, 
such as stra nge hyperons A° and S~, are expected to appear once their mass thresholds 
are reached (IBaldo et al.l 120001 ). The appearance of hyperonic degrees of freedom lowers the 
energy-per-particle in the stellar medium and causes more centrally condensed configura- 
tions with lower masses and radii. Additionally, hyperons help the condition for nucleonic 
direct Urea process to be satisfied at lower densities due to the increased proton fraction 
( iPrakash et al.lll992l ). and depending on their exact concentrations could pote ntially con- 



tribu te to the fast cooling of the star through hyperonic direct Urea processes (jPage et al. 



20061 ). These considerations would alter the balance between the curves in Fig. 9 and ulti- 
mately the results displaced in Fig. 8 in favor of the direct Urea process, i.e. smaller masses 
and higher frequencies would be necessary to close the fast cooling channel. This is due to 
the fact that the overall impact of (rapid) rotation on the neutron star structure is smaller 
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for more centrally condensed models resulting from "softer" EOSs (IFriedman et al.lll984j ). 
Therefore, for such models there is smaller deviation from properties and structure of static 
configurations. In addition, at even higher densities matter is expected to undergo a tran- 
sition to quark-gluon plasma (IWeberlll999l : iBaldo et al.ll2000l ). which favors a fast cooling 



throu gh enhanced nucleonic direct Urea and quark direct Urea processes (see e.g.. |Page et al. 
2006h . 



4. Summary 

We have studied properties of (rapidly) rotating neutron stars employing four nucleonic 
EOSs. Rapid rotation affects the neutron star structure significantly. It increases the max- 
imum possible mass up to ~ 17% and increases/decreases the equatorial/polar radius by 
several kilometers. Our findings, through the application of EOSs with constrained symme- 
try energy by recent nuclear terrestrial laboratory data, allowed us to constrain the mass of 
the neutron star in XTE J1739-285 to be between 1.7 and 2.1M0. 

Additionally, rotation reduces central density and proton fraction in the neutron star 
core, and depending on the exact stellar mass and rotational frequency could effectively 
close the fast cooling channel in millisecond pulsars. This circumstance may have important 
consequences for both the interpretation of cooling data and the thermal evolution modeling. 
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